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The thermodynamics of a scalar field with a quartic interaction is studied within the linear 8 
expansion (LDE) method. Using the imaginary-time formalism the free energy is evaluated up 
to second order in the LDE. The method generates nonperturbative results that are then used to 
obtain thermodynamic quantities like the pressure. The phase transition pattern of the model is 
fully studied, from the broken to the symmetry restored phase. The results are compared with 
those obtained with other nonperturbative methods and also with ordinary perturbation theory. 
The results coming from the two main optimization procedures used in conjunction with the LDE 
method, the Principle of Minimal Sensitivity (PMS) and the Fastest Apparent Convergence (FAC) 
are also compared with each other and studied in which cases they are applicable or not. The 
optimization procedures are applied directly to the free energy. 
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43 ' I. INTRODUCTION 

Oh| Phase transition phenomena in quantum field theories are typically of nonperturbative nature and thus naive 
perturbation theory based on an expansion in the coupling constant cannot be employed. This is clearly the case of 
phase changes at high temperatures, where perturbation theory becomes unreliable because powers of the coupling 
constant become surmounted by powers of the temperature Problems with perturbation theory also happen in 
phenomena occurring close to critical points, because large fluctuations can emerge in the system due to infrared 
divergences, thus requiring nonperturbative methods as well in their studies. This is the case of studies involving 
second order phase transitions and also in weak first order phase transitions (2j. Typical examples where these 
problems can manifest are in studies of symmetry changing phenomena in a hot and dense medium, a subject of 
interest in quantum chromodynamics (QCD) in the context of heavy-ion collision experiments, and also in studies of 
the early universe. Consequently, there is a great deal of interest in investigating thermal field theories describing 
matter under extreme conditions [1, 0, S ■ 

Familiar nonperturbative methods that have been used to study symmetry changing phenomena at finite tempera- 
tures are resummationlike techniques, such as the daisy and superdaisy schemes composite operator methods 0, 
and field propagator dressing methods (Tol . [llj • Other methods used include expansions in parameters not related to 
a coupling constant, like the 1/N expansion and the e-expansion [12J. In addition, there are numerical methods, the 
most notably ones are those based on lattice Monte Carlo simulations [HI . Each method has its own advantages and 
disadvantages. For instance, in numerical methods there may be issues related to numerical precision, lattice spacing, 
and lattice sizes. In addition, there is the notorious problem of simulating fermions on the lattice at finite chemical 
potentials [l3j . In any nonperturbative method based on an expansion in some parameter one has to face the problem 
of higher order terms becoming increasingly cumbersome, so stalling further analysis. This is usually the case when 
carrying out calculations beyond leading order in the 1/N expansion. Careless use of a nonperturbative method can 
also lead to problems like the lack of self-consistency or overcounting of effects. Known examples of such problems 
are the earlier resummation works dealing with daisy and superdaisy schemes, that at some point were giving wrong 
results, e.g. predicting a first order transition [l4| for the X<p 4 theory, an unexpected result since the model belongs 
to the universality class of the Ising model, which is second order. These methods also predicted a strong first order 
phase transition in the electroweak standard model, a result proved to be misleading |15j | . 
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Let us recall that the breakdown of perturbation theory at high temperatures and its poor convergence pro 



erties have been dealt with many different methods. Examples are the use of self-consistent approximations 
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hard-thermal- loop (HTL) resummation [17], [l8|], perturbative expansions in the coupling constant with resummation 
implemented with the use of a variational mass parameter, also known as screened perturbation theory (SPT) [T^.[2Cl|. 
and the use of two-particle irreducible (2PI) effective actions [2l|. The 2PI method, in particular, leads to a much 
better convergence of thermodynamic quantities (like the pressure) as compared to some of the other methods 22\ . 
Related to the 2PI method is the ^-derivable technique, which has been used to study the thermodynamics of scalar 
and gauge theories [23|, [H], HE, HE] • One difficulty with the 2PI effective actions is that the renormalization procedure 
is nontrivial (27j . In addition, there seems that the ^-derivable technique breaks down for a coupling beyond some 
value [H|. 

In general, it is desirable that any analytical nonperturbative method obey two basic requirements. First, it should 
be self-consistent, and second, it should produce useful results already at lowest orders without the need for going 
to higher orders. That is, it should produce results that quickly converge at some order where calculations are still 
feasible analytically or semianalytically. Though some of the cited methods may satisfy one, or to some extent both 
of these requirements, in the present paper we are particularly interested in the one known as the linear 8 expansion 
(LDE) 12911 . a nonperturbative method that has been used successfully in different contexts related to thermal field 
theories [30, HH, [32] and in many other theories - for a long, but far from complete list of references see Refs. [331 134]] . 
In the LDE, a linear interpolation on the original model Lagrangian density is performed in terms of a fictitious 
expansion parameter 8, which is used only for bookkeeping purposes and set at the end equal to one. The standard 
application of the LDE to a theory described by a Lagrangian density C starts with an interpolation defined by 

C^C 5 = (l-5)Co(v) + SC 

= Co (v) + 6 [C - C (77)] , (1.1) 

where £0 is the Lagrangian density of a solvable theory, which is modified by the introduction of an arbitrary mass 
parameter (or parameters) 77. The Lagrangian density C s interpolates between the solvable Co(r]) (when (5 = 0) and 
the original C (when 6 = 1). The procedure defined by Eq. (|1.1|) leads to modified Feynman vertices, that become 
multiplied by 5, and modified propagators, that now depend on 77. All quantities evaluated at any finite order in 
the LDE will then depend explicitly on 77, unless one could perform a calculation to all orders. Up to this stage the 
results remain strictly perturbative and very similar to the ones obtained via an ordinary perturbative calculation. 
It is through the freedom in fixing 77 that nonperturbative results can be generated in this method. Since 77 does 
not belong to the original theory, one may fix it requiring that a physical quantity calculated perturbatively to 
order-<5 fe , be evaluated at the value where it is less sensitive to this parameter. This criterion, known as the principle 
of minimal sensitivity (PMS), translates into the variational relation (35j 



drj 



= 0. (1.2) 

rj,d — 1 



The optimum value 77 which satisfies Eq. (jl.2p is a function of the original parameters of the theory. In particular, 77 is 
a nontrivial function of the couplings and because of this nonperturbative results are generated. Another optimization 
procedure used is known as the fastest apparent convergence (FAC) criterion [35|. It requires from the fc-th coefficient 
of the perturbative expansion 

k 

(1.3) 

that 



$(fc) _ $(*-!) 



6=1 



0, (1.4) 



which is just equivalent to taking the fc-th coefficient (at 8 = 1) in Eq. (|1.3[) equal to zero. 

One should note that it is not at all guaranteed that the condition in Eq. (| 1 . 2[) has a nontrivial solution. In cases 
where this may happen, the second criterion, Eq. (|1.4[) . may be more appropriate. One example where the condition 
given by Eq. (|1.2p fails to produce a nontrivial solution was in the problem studied by the authors in Ref . [36] , who 
applied the LDE to compute the effective potential in superspacc. There, the authors found that while the PMS 
condition was unable to give a nonperturbative solution to the effective potential, the FAC criterion worked perfectly 
well. Of course, in many situations both optimization criteria may work and in this case one may ask whether they lead 
to equivalent results. Previous studies indicated that this is indeed so, but a full comparison of results obtained with 
both optimization criteria is still lacking. Another issue associated with the LDE is its convergence. Rigorous LDE 
convergence proofs have been obtained for the problem of the quantum anharmonic oscillator, at zero temperature, 
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considered in Ref. [37|, while its partition function at finite temperatures was considered in [38|. For quantum field 
theories, Ref. [3^ | has proved convergence for a particular perturbative series in an asymptotically- free, renormalizable 
model at zero temperature. For a critical \<f) A 0{N) theory in three dimensions the issue of convergence was studied 
in [4(| employing both PMS and FAC optimization criteria. Finally, regarding the possible solutions that can emerge 
from the optimization criteria (PMS or FAC), we must use a definite approach in selecting the optimum root fj from 
either Eq. (|1.2|) or (|1.4|) . The problem of dealing with the many possible solutions for fj was treated in details in the 
first two papers cited in Ref. [4(|, where the convergence of the LDE was also studied in details. Typically, the higher 
the order in <5, the more solutions can appear. As shown in those references, all solutions at each given order in <5 
can be classified into families. The optimum value for r] is chosen as follows: The trivial solutions for fj, e.g. fj — 
and those that are not dependent on the coupling constant (and thus cannot lead to nonperturbative results) are not 
considered. In addition to these, at first order there is only one nontrivial solution (first family), consistent with all 
our approximations, (like the high-temperature approximation, used later in our calculations). This family is then 
followed in the next orders and used in all our calculations. As proved in earlier references with the LDE method, 
this is a consistent and unambiguous way for choosing the optimum value for rj. 

It is important to stress that in the method of the LDE the selection and evaluation of Feynman diagrams proceed 
in the same fashion as in ordinary perturbation theory, including the renormalization procedure [33l . l39l . |40(. The 
results obtained are free from infrared divergences, even at the critical point and in its neighborhood, thus making 
it a particularly suitable method to study phase transition phenomena in quantum field theories. It is important to 
recall here that there are similarities between the LDE and the SPT methods. In particular, the implementation of 
the latter can be put in a form similar to the LDE by means of a modified loop expansion [4l| . named optimized 
perturbation theory (OPT) in this reference. But there are also some major differences between these methods. For 
instance, in the LDE no assumption is made a priori for the parameter r\, while in the SPT/OPT it is assumed that 
such a mass term is already of some order in the coupling constant. The implication of this is that the order counting 
of loop expansion has to be readjusted accordingly. 

In the present paper we study the application of the LDE to the A0 4 theory. We will study the applicability of the 
PMS and FAC optimization criteria for the symmetric and broken phases of the theory, and compare results obtained 
with both methods. In addition, in the present work we choose to optimize the free energy, instead of the self-en ergy 
like in many other works employing the LDE, particularly Refs. [3l|,[32j]. There are several reasons for doing so (30l.[34j. 
but an important one is that in some situations it might happen that the optimization of the self-energy does not lead 
to nontrivial solutions, while optimization of the free energy with PMS or FAC are seen to lead to nontrivial solutions 
already at first order in S. The critical temperature T c , the pressure P, and the background dependent free energy F 
are obtained here in an explicit calculation up to order S 2 . Calculations at this order require a calculation of vacuum 
terms up to three loops. Since the thermodynamics of this model has been extensively studied before in the literature 
with a number of methods, our calculation here will be useful to benchmark the application of the LDE and its two 
main optimization procedures against those previous applications. In addition, we compare our results with those 
obtained with standard perturbation theory. Besides correctly reproducing the expected second order phase transition 
pattern for the model, our results at order S 2 are shown to be sufficient to obtain the thermodynamics of the model, in 
the sense that the results at 0(S 2 ) are not much different from the ones at 0(6). The results point towards a quickly 
convergent LDE, as already indicated in previous studies with different models under different conditions [39ll4fj|. 

This work is organized as follows. In the next section we introduce the interpolation procedure for the model. 
In Sec. Ill we compute the free energy in the symmetric and broken phases to 0(S 2 ). In Sec. IV we present the results 
obtained from the optimization procedures. The pressure is evaluated and contrasted order by order with the one 
obtained within perturbation theory. The critical temperature, the temperature dependent vacuum expectation value 
of the scalar field and the free energy are determined to 0(S 2 ). Our conclusions are presented in Sec. V. 



II. THE MODEL LAGRANGIAN DENSITY 



The interpolation defined by Eq. (jl.ip when applied for the standard A0 4 model gives 

C s = C (v)~S^ + S^ 2 + C s ct , (2.1) 

where 

Co(v) = \ (d^f - ^4> 2 - ^ 2 , (2.2) 

and C,\ t is the part of the Lagrangian density carrying the renormalization terms needed to render the model finite. 
Details about this renormalization procedure in the LDE and the explicit form for C s ct are given e.g. in Ref. (3lj for the 
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case of background field dependent contributions (broken symmetry phase), while the field independent contributions 
(symmetric phase) were given in Ref. [20| within the context of the SPT, so we will not repeat those same renormal- 
ization details here. One should also note that the only "new" terms introduced by the 5-expansion interpolation are 
quadratic terms and so the renormalizability of the original theory is not changed. This means that the renormaliza- 
tion of the theory can be carried out in an analogous way as in ordinary perturbation theory [3~fl ] . Specifically, the 
interpolation procedure given by Eq. (|2.1j) introduces a new (quadratic) interaction term, with Feynman rule iSrj 2 . In 
addition to this modification, the original bare propagator, 

S(k) = i (k 2 - m 2 + ie)' 1 , (2.3) 

now becomes 

S s (k) = i (k 2 - ml - rf + ie)' 1 , (2.4) 

while the original quartic vertex is changed from — iX to — iSX. 

In the next section we will show the results for the finite temperature free energy density using the interpolated 
model with the LDE at 0(6 2 ). We will consider the cases of m 2 = \m^\ 2 and to 2 , = — |mo| 2 in Eq. (|2.2[) . corresponding 
to the symmetric and broken phases, respectively. 



III. THE FINITE TEMPERATURE FREE ENERGY IN THE LDE TO 0(S 2 ) 

We perform the standard derivation of the free energy [i^] up to 0(S 2 ). With the constant field introduced through 
the usual shift of the scalar field, 4> — ► 4> + ip, the Lagrangian density is rewritten as 

£[4>(x),Lp]=£2[<t>(x),<p}+/: I [4>(x),<p} , (3.1) 

where £2 is the part of the Lagrangian quadratic in the fields, 

1 ,„ ,,2 1 



while Ci is 



where in Eq. (|3.2p il 2 is given by 



L 2 ^{x),ip]^-{d^Y , (3.2) 



£zfe(s),d = -y^ 3 -^ 4 , (3.3) 



n 2 = ±m 2 + S -± ip 2 + (l-5) V 2 . (3.4) 

Note that in all loop contributions the propagators will carry a mass term as given by Eq. (|3.4[) . These terms are then 
expanded in 5 to the desired order, thus generating the insertions of rj 2 that appear as a consequence of the quadratic 
vertex introduced in Eq. (|2.1[) . 
The free energy is 

F[tp] = F (ip) + F^ loop (p) 

+ ^iln(^exp l^i J d*x d[(l)(x),Lp]^ \ , (3.5) 

where Fo(ip) is the tree-level classical potential and Fi-\ oop (ip) is the one-loop contribution to the free energy (V is 
the space volume) given by 

Fi-ioop(^) = y i In fdcf> jJ^lrM*)] . (3 . 6) 
Higher loops arc given by the last term in Eq. (|3.5[) . with the average (• • • ) meaning 

fV<j> (■■■) e ifd i xC 2 l<p(x)M 



§V(j> e ifd 4 xC 2 [4>(x),<P] 



(3.7) 
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As said above, the scalar field propagators in the diagrams are obtained from Ci\(j>{x), if], and the vertices are 
determined from Cj[(j){x), tp], with both as given at the end of Sec. II. 

Our calculations are performed, as usual, in the imaginary-time formalism (42j. Thus, the scalar boson field has 
Euclidean four- momentum P — (u n , p), with P 2 = ui^ + p 2 , where oj n are the discrete Matsubara bosonic frequencies 
Lj n = 2im/f3, with n = 0,±1,±2, • • •, and j3 = l/T. Loop diagrams involve sums over the Matsubara frequencies 
and integrals over the space momentum p. All space momentum integrals are performed in arbitrary dimension 
d = 3 — 2e and renormalization is performed in the modified minimal subtraction scheme (MS). The measure used in 
the sum-integrals is then defined as 



e> 2 

47T 



(27T) 



3-2e ' 



(3.8) 



where fi is an arbitrary momentum scale in dimensional regularization. The factor ( 6 4 ^ ^ is introduced so that, 

after minimal subtraction of the poles in e due to ultraviolet divergences, [i coincides with the renormalization scale 
in the MS scheme. 

From Eq. (|3.5p . the free energy is expressed up to O (<5 2 ) by expanding all appropriate terms in 5. Considering 
the vacuum contributions to the free energy, this means that terms up to three-loops must be included. All bare 
(unrenormalized) contributions are shown in Fig. [1] 



(a) 









(b) 





Figure 1: Diagrams contributing to the free energy up to 0(S 2 ), given by (a) vacuum diagrams and (b) background field 
(external legs). The black dots indicate a Srj 1 insertion. 



The renormalization procedure for the symmetric phase was performed in detail in Refs. [2fJ,|31(. The counterterms 
for the vacuum diagrams are given in Ref. (2pj |. while those for the field dependent diagrams are given in Ref. [3lj |. 
We also note that the di verg ences in the broken phase can be removed by the same counterterms determined for 
the symmetric phase [4l|, |43[ HH, so the renormalization for the broken phase does not require extra effort. The 
renormalization proceeds just as in standard perturbation theory and as shown in detail in Ref. (3lj |. only tempera- 
ture independent counterterms are required and the temperature dependent divergent terms cancel out exactly. All 
diagrams of counterterms contributing to F[ip] up to 0(6 2 ) are shown in Fig. O 

(a) ( (^^) ^C^^ 6 




(b) 





Figure 2: Diagrams representing the counterterms for the free energy up to 0(\ 2 ): (a) vacuum contribution, (b) background 
field contribution. As in Fig. 1, the black dot indicates a Srf insertion. The circle-cross denotes either insertion of a mass 
counterterm or of a vertex counterterm. 



The circle-cross in Fig. [5] denotes either a mass counterterm vertex Am 2 , or a vertex counterterm AA, given 
respectively by [3~fl | 



Am 2 
AA = 



<5 32k^ 2 + (1 -^ 2 )]-^ A " 
32tt 2 £ 



(32tt 2 ) 2 



-2 1 

T + - 

r e 



(m 2 + rj 2 ) 



(3.9) 
(3.10) 
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The final expression for the renormalized free energy F[ip], including all terms shown in Figs. Q] and [2] becomes 



-^l^] — ^vacuum ~t~ Ftp 



(3.11) 



where F v , 



denotes the vacuum contributions, 



1 



vacuum 
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-Jo{f3M)T 4 
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+ 1 M* - J\ (PM)T< 
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A 



4(4tt) 4 
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48 (4vr) 



ln ( T77 ) + J 2 



Inf^J+lJ^-Jx^)^ 



51n 3 f^f )+171n 2 



f (l) + C* + 3 In 



M 2 



+ 1 J 2 (/3M) 



41, / M 2 
— ln 

12 In 



M 4 



23 



21M 2 



23 
12^2 



28 ln 



12 - 7T 2 - 4Ci + 6 ln 



1 J 2 



Ji (f3M)n 2 T 2 



+ 



3 31n 



/x 2 



+ 4 J 2 (/37W) + 3Jf (pM) J 2 (#M) + 6K 2 + 4K 3 



T 



(3.12) 



and F^ denotes the background field dependent contributions, 



Fin — F 



f SXM' 2 
A 2 M 2 



l0g 1 If 



, At? 2 
32tt 2 



I ^- ) ;.v() 



2(32tt 2 )' 



In 
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M 2 



? \ \ 2 o 



, 3A 2 X 2 
2(32tt 2 ) 2 



In ( ^ ] - I j + I 



7T 

IT 



1024tt 4 



M * ( 1 + — ) + 4 - A/[2 ln (ja) I 1 + J 2 (z 3 ^)] + J 2 ^ 



8M 2 In 2 (ii) - 4 ln (ii) J x (/?M) T 2 - J 2 (/3A4) J x (/?M) T s 

p A 2 r 2 

24 (4tt) : 



.M 2 

In I I + : <"' !(, - >r > 



2 A 2 M 2 
256tt 4 



24 V H 



2 In 2 ( — I + 1.164032 



ip- 2 3A- 



2 5 32 tt 



(3.13) 



with i^o given by 



F = I [±ml + (1 - *)r? 2 ] ^ 2 + « 



(3.14) 



In Eqs. (j3~T2f and (J3TT3J) , .M 2 = im 2 , + ry 2 , and the constant terms appearing in Eq. (|3 . 1 2[) are defined as follows: 
^"(1) = -2C(3), where ((x) is the zeta function, C 2 ~ 39.429 and C 3 ~ -9.8424, while K 2 and K 3 are three- 
dimensional integrals that can be evaluated numerically [ia ]. In the high-temperature limit, M/T <C 1, they are 
given by 20] 



32tt 4 



[ln (#M) - 0.04597] - 372.65/37W [ln (/3M) + 1.4658] 



(3.15) 
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and 



K 3 ~ 453.51 + 1600/iM [In (/37W) + 1.3045] . (3.16) 
In Eqs. (|3.12p and (|3.13p . we have also denned the temperature dependent integrals J n (n — 0, 1, 2) as follows, 

c 4-2„ j 



J n (a)= [2 > 



dx- 



r (| - n) 7 %/z 2 + a 2 e v / ^+? 
which can be expressed as a series expansion as follows [l|, H2, 0, [47| 

8tt 



1 



J (a) 



-a d + a 4 In 



(Ks)^"!) 



128 £ 



(-l) n (2n- l)!!C(2n + l) a ( 2 ™+ 4 ) 



32 (n + 2)!2"+ 1 (27r)' : 



4tt 2 o 16 4 
3 45 



(3.17) 



(3.18) 



Ji(a) = -47ra - 2a 2 
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and 



J 2 (o) = 



2tt „, / a \ 
— + 21n — 
a V47T/ 



2 7 



E 



(-l) n (2n- l)!!C(2n+l)a 2 
n!2"+! (27r) 2 ™ 



(3.19) 



(3.20) 



Equations (|3.18| - (|3.20p are all convergent in the high-temperature limit as can be easily checked by considering a few 
terms in the sums in these equations. 

We should note that when optimizing the free energy, since Jo, J\ and J2 are dependent on 77, it is important 
to check the stability of the results when truncating the sums in Eqs. (j3. 18|) - ()3.20|) . This is particularly critical for 
parameter values such that M./T is not much smaller than 1, a situation that requires a fairly large number of terms 
in the sums. In all results shown in the next section we have used enough terms in Eqs. (|3.18[) - (|3.20p so to obtain 
stable results for all parameter and temperature values used. 



IV. OPTIMIZATION AND NUMERICAL RESULTS 



We now turn to the application of the optimization procedures in the LDE and show the results obtained by 
implementing the PMS, Eq. (jl.2|) . and FAC, Eq. (|1.4p . As we explained in the introduction, the optimization criteria 
are applied directly to the free energy. The results obtained with each optimization criterion are contrasted with 
each other and with those available from other methods. This will then allow us to gauge the performance of each 
optimization procedure regarding both reliability and convergence. 



A. Symmetric Phase 

We initially restrict our calculations to the symmetric phase (with positive mass term in the classical potential) and 
evaluate the pressure, P = —F. In Fig. [3] we show our results for the pressure using the usual perturbation theory in 
A up to C(A 2 ) and where we have restricted to the case of high temperatures (M /T<1). In this figure the behavior 
of the pressure is shown as a function of the renormalized coupling constant, and T Q is a reference energy scale 
chosen as T = mfl(T ), where is the renormalized mass. This is similar as done in Ref. [22| using the 2PI method. 
Note that in Ref. [22| the authors define the quartic coupling differently from us. In their case, g 2 = A/24, and their 
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results are plotted as a function of the renormalized coupling g^. We choose here the same scale as in Ref. [22[ so to 
facilitate the comparison with their results for the pressure. It is clear in Fig. [3] the typical alternating behavior of 
the perturbative calculation, which indicates its very poor convergence. 



1.02 
1.01 
1 

n^O.99 
0.98 
0.97 

°' 96 0.2 0.4 0.6 0.8 1 

(X R 124 ) 1 ' 2 

Figure 3: The pressure in the symmetric phase using perturbation theory in A up to 0(A 2 ). The parameters used are T = 20 To 
for the temperature and fi — To for the renormalization scale. 

Next, we use the result for the free energy evaluated up to 0(S 2 ), given by Eq. (|3.12[) . Note that in the symmetric 
phase the pressure depends only on vacuum terms, since the free energy is evaluated at the vacuum expectation value 
for the field, <p = 0. By optimizing the free energy using the PMS criterion, Eq. (|1.2|) . we determine the root fj, which 
is then substituted back into the expression for the free energy, with the criterion used for choosing the optimum 
root as discussed below Eq. (|1.4[) . This naturally brings nonlinear A contributions and generates nonperturbative 
results. The pressure obtained in this case is shown in Fig. [4j where we show the results obtained up to orders S and 
5 2 . In the same figure we also show the perturbative results of Fig. [3] for comparison. It becomes evident here the 
convergence of the results with the LDE-PMS, with both 0(5) and 0(S 2 ) results not differing too much, in contrast 
to the perturbative (in A) results. 



i 1 i 

— perturbative X. 1 

2 

perturbative \~ 




1.02 r 



1.01 



perturbative X 
perturbative X' 
PMS 8 1 
PMS 5 2 




0.6 
(V24) 1/2 



Figure 4: The pressure in the symmetric phase to orders S and 8 2 using the PMS optimization criterion 
results of Fig. [3] are also shown. The parameters used are the same as those in Fig. [3] 



the perturbative 



In Fig. [5] we show again the results for the pressure, but now using the FAC optimization criterion, Eq. (|1.4[) . We 
once again see the excellent convergence for the pressure when contrasted to the perturbative results. 

In Fig. [6] we plot side by side our results for the pressure at order 5 2 using the PMS and FAC optimization criteria. 
It is seen as an excellent agreement between the two optimization criteria and it shows the equivalence of these two 
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Figure 5: The pressure in the symmetric phase to orders S and 5 2 using the FAC optimization criterion, 
optimization procedures. 
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Figure 6: The pressure for the symmetric phase at order 5 2 using the PMS and FAC optimization criteria. 

A side by side comparison of the order 6 2 result for the pressure (from either the PMS or FAC) with the 2PI 
two-loop result of Ref. [22| (second panel of theirs Fig. 1) shows an excellent agreement between the two results. Since 
operationally the LDE is much simpler to be implemented than the 2PI calculation and also when compared with 
other methods, like those based on the renormalization group and Schwinger-Dyson equations, this may be a great 
advantage of the LDE. Many previous applications of the LDE to a large variety of problems (cited previously) also 
confirm the strength of the method. Its strength comes basically from the fact that its implementation is similar to 
that of standard perturbation theory. The important and fundamental difference with standard perturbation theory 
resides in the optimization procedure that fixes an initial, a priori, arbitrary parameter, r\. It is then interesting to 
investigate what kind of role the optimum r\ represents in the LDE after optimization. This is partially clarified in 
the plot shown in Fig. [7J where we show the optimum rj as a function of the renormalized coupling constant. It shows 
that by increasing the order in 5, fj becomes closer and closer to the thermal mass my, here computed at one- loop 
order for simplicity. In general, we can extrapolate this expectation and say that the expected optimum value for r\ 
should be the (quantum and) thermal mass (quadratic in the field) corrections, as would be derived from a true gap 
equation. This is in fact confirmed by the many applications of the LDE to the Gross-Neveu model [30], in which case 
exact results are known (in the large- A" approximation) and can then be readily compared with the results obtained 
from the LDE method applied to that model. 
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Figure 7: The behavior of the optimum parameter rf with respect to the renormalized coupling constant and evaluated at 
order 5 and 5 2 using the FAC optimization criterion. The parameters used are the same as in the previous figures. 



B. Broken Phase 



Let us now turn to the symmetry broken case (with negative square mass term in the classical potential). For 
this case we found that only the FAC optimization criterion leads to nontrivial solutions for r\. The FAC criterion is 
applied to the free energy and the resulting nonlinear equation is solved simultaneously with the equation giving the 
minimum condition for the field (thermal) expectation value, v{T) 1 given by the condition 



dF[ip] 



d<p 



= 0. (4.1) 



As it is well known, the phase transition in the pure scalar theory is second order [T^], as required by universality 
reasons. Our results for the free energy using the FAC criterion indicate a second order phase transition. This is 
shown in Fig. [Sj where the free energy for A = 0.1 is plotted for different temperature values. The critical temperature 
obtained here is T c j /i ~ 15.49, consistent with the perturbative prediction |42| and other nonperturbative calculations 
[To| . Another quantity that indicates that the transition is a continuous one is the temperature evolution of the 
minimum of free energy, v{T). This is shown in Fig. [9] for A = 0.1 and A = 1.0. 
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Figure 8: The nonperturbative free energy for A = 0.1 evaluated at order 5 2 and using the FAC optimization criterion, for three 
different temperatures: T < T c , T = T c and T > T c , where T c = 15.49 (in units of the renormalization scale fi). Here we have 
set [i = mo. 
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Figure 9: Temperature dependence of the minimum of the free energy v(T) at order S 2 obtained with the FAC optimization 
criterion. Here u(0) is the tree-level minimum of the bare free energy. 



Finally, in Fig. [10] we show the temperature dependence of the thermal mass, mx, as derived from the free energy, 



d 2 F[<p] 



d 2 tp 



1,3=0 



(4.2) 



We once again can notice a continuous and smooth transition. We note that one can determine the critical temperature 
by looking at which value of T m 2 ^ changes sign and check whether this gives the same result for T c as obtained from 
v(T c ) = (as in Fig.©. 




Figure 10: Temperature dependence of the nonperturbative thermal mass at order 8 2 evaluated with the FAC criterion for two 
values of the coupling constant A. We use in this plot fi — mo- 



V. CONCLUSIONS 



One of the motivations for using the LDE to study the thermodynamics of the scalar field theory at high tempera- 
tures, as done in this work, was its ease of implementation and renormalization, which is no different from those of the 
standard perturbation theory. One recalls that similar studies in the context of the 2PI and related methods typically 
face difficulties in the renormalization procedure, making their applicability a nontrivial task. In addition, the LDE 
method, differently from other methods, like SPT (or OPT - optimized perturbation theory), makes no assumption 
on the introduced mass parameter rj, thus we do not need to adjust the order of the loop expansion accordingly. 
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By using the LDE, we have studied the thermodynamics of the \<p A scalar field theory in the symmetric and broken 
phases. The LDE is used with two popular optimization procedures, known as PMS and FAC. There are two major 
differences with the work we have done here and previous ones, like e.g. in Refs. [3TI. I32I]. First, while in general the 
PMS procedure is the favorite optimization criterion in the literature related to the LDE and similar methods, we 
here have shown that the FAC procedure leads to numerically indistinguishable results from the ones obtained with 
the PMS. In addition, while there may be cases where the PMS procedure leads to trivial results only, the same may 
not be the case for the FAC (here we have shown this to be the case in the broken phase). In this sense they can 
be used in a complementary way, when PMS fails, one can try FAC, or vice- versa. Secondly, unlike in Refs. [3lL l32|. 
where the quantity optimized is the self-energy, here we choose to optimize the free energy. One advantage of this is 
that, while there is no solution for the LDE at first order when optimizing the self-e nerg y, we do find solutions when 
optimizing the free energy already at first order in 5. Furthermore, as shown in Ref. [301 ] . the optimization of the free 
energy can be shown to immediately lead to the solution of the gap equation (here verified numerically through the 
results for the optimum 77), while in optimizing the self-energy further constraints must be employed, as for example 
renormalization group equations. In the numerical studies performed in the present work, we have shown that the 
optimum 77 carries both temperature and coupling constant contributions. Thus, the LDE with optimization of the 
free energy implements automatically a nonperturbative resummation of the thermal corrections, in conformity with 
analytical results produced when this method was used to study the Gross-Neveu model in Ref. [30j , from which exact 
solutions are available and a close comparison with the LDE results is possible. 

By studying the behavior of the pressure and contrasting the results obtained with perturbation theory and the 
2PI method, we have shown that the LDE leads to convergent results already at lowest order in the LDE expansion 
parameter 5, with the first and second order results changing only slightly and producing results consistent with 
the 2PI nonperturbative method. In addition, as already mentioned above, we have shown that both optimization 
procedures, FAC and PMS, lead to equivalent results for the pressure. 

Another important result of our work is that the LDE is shown to be adequate for studying phase transitions at high 
temperatures. In particular, when applied to the phase transition in the A</> 4 , the LDE predicts the correct order of the 
phase transition, which is second order, in agreement with general results of statistical mechanics. Besides this, since 
the LDE method automatically produces an infrared cutoff, the results are shown to be valid and applicable below, 
above, and at the critical temperature T c , showing that the LDE circumvents the usual problem seen in perturbation 
theory, namely, the appearance of infrared divergences close to critical points. 
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